Full counting statistics of electron transfer between superconductors 
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We present an extension of the Keldysh-Green's function method, which allows to calculate the 
full distribution of transmitted particles through a mesoscopic superconductor. The method is 
applied to the statistics of supercurrent in short contacts. If the current is carried by Andreev 
bound states the distribution corresponds to switching between long trains of electrons going in 
opposite directions. For weak (gapless) superconductors or tunnel junctions we find that at low 
temperatures the distribution has negative "probabilities" . Accounting for the quantum mechanical 
nature of the measuring device shows that these negative values can indeed be measured. 
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Coherent charge transfer between superconductors (S), 
supercurrent, is essentially a quantum- mechanical pro- 
cess. Although superconducting junctions are commonly 
used, the statistical properties of the charge transfer in- 
volved into the supercurrent are not yet completely un- 
derstood. In view of recent attempts to use the coher- 
ence of superconductors to build quantum bitsQ], it is 
necessary to reveal the basic limitations on this coher- 
ence (if there are any) . Additionally, the problem of the 
statistics of transferred charge in a quantum process is of 
fundamental interest. It is related to the understanding 
of the measurement process and the interpretation of its 
outcome. 

Recently the current noise exhibited in SNS junc- 
tions, where N is a diffusive normal metal, was addressed 
experimentally. Q The experimental results show a giant 
excess noise in the low temperature and voltage regime 
in those samples, in which at the same temperature a 
coherent coupling through the normal metal was mea- 
sured. This is in accordance with theoretical predictions 
for the shot noise in short contacts. || This may hint 
to the importance of an understanding of the statistical 
properties of the supercurrent in such junctions. The 
equilibrium noise properties have been studied in jm and 
p). Further experimental progress in the fabrication of 
controllable single-channel junctions is to be expected in 
the near future. This will shed more light on the funda- 
mental statistical properties of charge transfer between 
superconductors. 

We will make use of the so called full counting statis- 
tics (FCS), originally introduced to calculate the distri- 
bution of transmitted charge through a contact between 
normal metals This theory allows to find the cumu- 
lant generating function (CGF) S(x), from which the 
distribution of transmitted charge follows via P(N) = 
J dxexp(—S(x) — iN\)- It is tempting (and has been 
done so far) to interpret P{N) as the probability that N 
charges have been transferred through the contact dur- 
ing the time of observation. We will show below that 
this interpretation is strict only for normal constrictions. 



For superconducting constrictions the distribution also 
depends on the phase difference <p. It turns out that 
P(N, 4>) can also take negative values, which hampers 
such interpretation. This is related to the fact that the 
phase and number of charges transferred can be regarded 
as canonically conjugated variables. Still P(N, 4>) pro- 
vides a complete description of all charge transfer pro- 
cesses and can be extracted from the results of measure- 
ments. 

The most powerful and general method of calculat- 
ing transport properties of mesoscopic conductors is the 
nonequilibrium Green's function approach (see 0). It 
was shown in Q that this approach can be generalized 
to access FCS. In this Letter we extend the approach to 
superconductors. This allows us to obtain the FCS of an 
arbitrary mesoscopic conductor at all temperatures and 
voltages. The CGF is derived for a contact which is fully 
characterized by an ensemble of transmission eigenval- 
ues {T n }. We evaluate the FCS of supercurrent in two 
generic cases. First, we find the distribution of trans- 
mitted charge of a single channel contact between two 
gapped superconductors. Here the current is carried by 
phase dependent Andreev bound states and, as shown 
by our analysis, conforms with the switching picture.^] 
The two bound states carrying current in opposite di- 
rections are alternately occupied and charges are trans- 
fered in 'long trains', which reflects the coherent nature 
of the supercurrent. In the second case we calculate the 
CGF of a contact between two weak superconductors. 
The resulting CGF corresponds to the tunnel limit for 
gapped superconductors and can be related to the ef- 
fective Keldysh action of a Josephson junction discussed 
in detail in Ref. ||. The standard interpretation^ of 
the CGF leads in the low temperature regime to nega- 
tive "probabilities" P(N,<fi). Negative values of P{ N, <f>) 
occur because of an attempt to interpret the quantum 
mechanical phenomenon of supercurrent with classical 
means. If we account for the quantum mechanical na- 
ture of the measuring device, we can resolve the paradox 
and specify how P(N, (f>) can be measured. 
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To be specific let us now introduce our model system, 
which is depicted in Fig. [j]. A mesoscopic conductor is 
placed between two reservoirs. The counting field x is 
introduced on an arbitrary cross section in one of the 
reservoirs and couples to the operator of current through 
that cross section. It follows from the definition of the 
cumulants that the CGF can be found from 
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Here T(T) denotes the (anti)time ordering operator. / 
denotes the current operator J d 3 x4f ji f 3 (/m)^ ji \/F(x), 
where 4 r is the usual Nambu spinor field operator and 
T3 is a matrix in Nambu space. Vi 7, is chosen such, that 
the spatial integration is restricted to the cross section 
and yields the total current. The counting field param- 
eterized in this way can now be incorporated into the 
boundary condition imposed by the reservoir onto the 
mesoscopic conductor. || That is, the reservoir Green's 
function effectively takes the form 



(2) 



Here Gi(</>) is the reservoir Green's function at super- 
conducting phase 4> i n the absence of the counting field 
and tk = &3T3 a matrix in Keldysh(")-Nambu(~) space. 
Now the counting field is included into the boundary con- 
dition for the Keldysh-Nambu matrix Green's functions 
provided by the left reservoir. Inside the system of inter- 
est the transport properties are described by quasiclas- 
sical Eilenberger equations [Q, applicable if the system 
size exceeds the Fermi wave length. It is important that 
G(x, 0) still obeys the quasiclassical normalization con- 
dition G 2 (x» = 1. 
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FIG. 1: Sketch of the system. Two reservoirs (1,2) are con- 
nected to a mesoscopic conductor M. The counting field \ ls 
chosen nonzero on the cross section C in reservoir 1. 



In certain cases the action S(x, 4>) can be found quite 
generally. One example is a constriction shorter than the 
coherence length, which is fully characterized by a set of 
transmission eigenvalues {T n }. The counting field ma- 
nipulates the matrix structure of the Green's functions 
in Keldysh-Nambu space. To find the transport proper- 
ties one should therefore use expressions, which respect 



the full matrix structure. It was noted in [|T(|, that a 
convenient way to do this is to use a "matrix current" , 
which is conserved in short contacts. The matrix current 
is formed with the current operator and the correspond- 
ing matrix elements of the Green's functions. Physical 
currents are related to certain components of the matrix 
current. For our purpose here, we need the matrix cur- 
rent in a short contact. The matrix current was derived 
in |l(J in the absence of the counting field. The counting 
field does not change the matrix structure of that result. 
So we can use it just by including the x _ dependence of 
the Green's functions and write 
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+ T n ({G 1 ( X ,0),G 2 }-2) 



(3) 

The action can then be found from the relation 
(d/dx)S(x,<t>) = -it Tr(f K i(x, Using the fact that 
[A, {A, G2}] = for all matrices with A 2 = 1, it is easy to 
verify that under the trace in (||) (d/dx){Gi(xi 4>)i ^2} = 
ifx [Gi(x, (j>), G2] . We can therefore integrate Eq. (||) 
with respect to x and obtain 

S(X,4>) = ^Y. I ^Trln[4 + T„({G 1 (x»,G 2 }-2)] 

n 

(4) 

Eq. (|4|) is very general. It contains the statistical prop- 
erties of all types of superconducting constrictions. For 
instance the FCS of an SN-contact can easily be ob- 
tained from 

In the rest of the paper we will study equilibrium noise 
and statistics of systems with two superconducting con- 
tacts. We will distinguish two generic cases. The first 
will be a single-mode contact of arbitrary transparency 
between two fully gapped superconductors. In the sec- 
ond case we treat a contact between two weak supercon- 
ductors, or, which is equivalent, a tunnel contact. The 
channel summation in the action (|j) is then a trivial sum- 
mation over transparencies. In the following derivation 
we limit the discussion to a single channel of transmission 
T\ and identical reservoirs at equilibrium. To be specific, 
we consider the Green's functions of the reservoirs 
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Here R(A)(E) are retarded and advanced Green's func- 
tions of the banks and h(E) — tanh(i?/2T) accounts for 
the equilibrium distribution at a temperature T . The 
phase difference </> is introduced by setting G\{<p) — 
exp(i0f 3 /2)Gs exp(— i0f 3 /2) and G2 = Gs. Advanced 
and retarded functions in (||) possess the structure 
-R(^4) = gK,AT% + /r.ati fulfilling the normalization con- 
dition f 2 + g 2 = 1. They depend on energy and the 
superconducting order parameter A. Their precise forms 
will be defined below. 
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The trace in the action can be evaluated and we obtain 
the main result of this paper 



S{X,4>) = ^ J dE\n 



1+E 



M<t>) 
Q{<t>) 



(e mx - 1) 



(6) 

Introducing q = (1 - 3r<?a)(1 - h 2 ) + /r/a(1 + h 2 ) the 
coefficients may be written as 
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(8) 



.(9) 



The interpretation of the different terms is analogous to 
that given in |12] ], A coefficient A± n is related to events 
in which a charge n is transfered to the right(left). The 
presence of terms, which describe charge transfers of 2e, 
is a consequence of superconducting correlations. The 
interpretation of these terms as probabilities stems from 
the comparison with the case of binomial statistics (see 
||). As we will discuss below this interpretation only 
makes sense for normal metals. 

Considering only the phase independent terms of (^])- 
(^) demonstrates an interesting feature. These terms can 
be factorizcd into 



T iq 



1 



(10) 



Due to the logarithm in Eq. (||), the exponent of 2 can 
be written as prefactor to the CGF. Therefore, the CGF 
describes two statistically independent probabilistic pro- 
cesses, which we can identify with electron and hole 
transfers. The phase dependent terms inhibit this fac- 
torization and lead to correlations between electrons and 
holes, as expected in superconductors. The denomina- 
tor Q, common to all coefficients, has roots for energies, 
at which Andreev bound states exist. Consequently the 
statistical properties are dominated by the charge trans- 
fer through Andreev bound states. The corresponding 
"probabilities" can be very large, in particular, larger 
than 1. Thus, we can not interpret the coefficients ([?])- 
@ as probabilities anymore. To find the statistics of the 
charge transfer, we have to specify the system further. 

Gapped superconductors. If the two leads are gaped 
like BCS superconductors the spectral function are given 
by /r,a = iA/((E± iS) 2 - A 2 ) 1 / 2 and g KA follows from 
normalization. Here S is an broadening parameter, which 
accounts for the finite lifetime of the Andreev bound 
states due to e. g. phonon scattering. The supercurrent 
is solely carried by Andreev bound states with energies 



±A(1 - Tisin 2 0/2) 1 / 2 = ±£ B (</>)• The importance of 
these bound states can be seen from the coefficient Q (|j) . 
It may become zero and will thus produce singularities in 
the action. |l3j The broadening S shifts the singularities 
of Q into the complex plane and allows an expansion of 
the coefficients in Eq. (0)-@ close to that energy. Per- 
forming the energy integration the action results in 



S(X, 



= -2*, 



tfyJl-IfttW/W-ixIiW/S, (11) 



where I\{4>) = A 2 Ti s\a{4>) /2E^{4>) is the super- 
current carried by one bound state and Ii(<p) = 
Ii(4>) tanh(£ , B/2T) is the average current through the 
contact. In deriving ( [TT| ) we have also assumed that 
X« 1. This corresponds to a restriction to "long trains" 
of electrons transfered, and the discreteness of the elec- 
tron transfer plays no role here. Fast switching events 
become less probable at low temperatures and are ne- 
glected here. In the saddle point approximation at low 
temperatures 7 = 1/ cosh(i?B/2T) « 1 we find for the 
current distribution 

Ptj 0) ~ ^^{i^pw-iwV^) r 12 ) 
7 

for \j((f>)\ < 1 and zero otherwise. Here, we have ex- 
pressed the transfered charge in terms of the current nor- 
malized to the zero temperature supercurrent: — 
I/Ii(4>). The current is related to the particle number 
by N — Hq. At zero temperature Eq. ( |l2] ) approaches 
P(j,</>) — > S(j — 1), which follows from a direct calcula- 
tion. Thus, at zero temperature the current is noiseless 
and the distribution (|l2|) at finite temperature confirms 
the picture of switching between Andreev states which 
carry current in opposite directions, suggested in Ref. [g]. 

Let us finally comment on the limits under which the 
previous result is valid. In the energy integration it was 
assumed that the bound states are well defined. For small 
transmission the distance of the bound state to the gap 
edge is w J\A. Thus, to have well defined bound states 
we have to require 5 < TiA. Similarly for a high trans- 
missive contact and a small the phase difference we re- 
quire (f) ~ I/I c > S/A. The statistics beyond these limits 
is similar to what is discussed in the following. 

Tunnel junction/ gapless superconductors. Let us now 
consider the supercurrent statistics between two weak su- 
perconductors, where the Green's functions can be ex- 
panded in A for all energies. One can see that this is 
equivalent to the tunneling limit ({T„} <C 1) of Eq. (Q). 
We also return the many channel situation here. Ex- 
panding the action (Q) to lowest order and using that 
the counting rotation can be written as exp(i%fif/2) = 
1 [e l x/2(! _ fjf ) + e -ix/2(j + f K )] we find 



S = -t [il s ((f>) sin X + P s {4>) (cosx - 1)] 



(13) 



In short, the full statistics are expressed in terms of su- 
percurrent is(<^) and noise P s {4>)- In equilibrium using 
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-Re J dETr {f 3 [iZifa), iia] } h (14) 
-Re J dETr {f 3 Ai(0)f 3j R 2 } (1 - fr 2 )(15) 



Here G = (l/V) T n is the normal state conductance of 
the contact. The equivalence of this result to the limit 
of gapless superconductors follows from an expansion of 
(@)-(|) to orders f 2 . 

Eq. (|l5|) shows that P s vanishes at zero temperature, 
since h(T = 0) = ±1, whereas I s vanishes at T c . There- 
fore, there is some crossover temperature below which 
P s < I s - In this limit the action possesses no saddle 
point anymore and by expansion in powers of exp(ix) it 
follows, that P{N, 0) becomes negative. Obviously this 
questions the direct interpretation of P(N, (f) as a prob- 
ability. Thus, we are forced to have a closer look on what 
P(N, 4>) actually is. 

To clarify this issue, we make use of the recent results 
presented in |OL where it was shown that the interpre- 
tation of P(N, (f>) in intimately related to the way the 
measurement is performed. We assume a simple model 
of a measuring device: a capacitor of infinite capacitance 
that stores the charge passed through the constriction, 
i. e. the charge operator q is related to the current op- 
erator through the constriction by q= I. The quantum 
mechanical treatment of this device involves its density 
matrix p{q 1 q'). 

In ]l4| the relation between initial and final density ma- 
trices of the device was obtained. This can be expressed 
in terms of the density matrix in Wigner representation, 
p(x,q), x being the canonical conjugate of q. It reads 

Pf{x,q)=Y J P{N^~x)p l {x,q-N), (16) 

N 

so that P(N, cj>) fully characterizes the quantum mechan- 
ical behavior of the capacitor. For a normal constric- 
tion P(N, (f>) does not depend on <j). In this case we can 
rewrite Eq. ([l6|) directly in terms of charge distributions 
U(q) = J dxp(x,p), 



n f (q)=J2P(NWq-N). 



(17) 



N 



Therefore, P(N) can be interpreted as classical proba- 
bility. For a superconducting constriction quantum me- 
chanics is essential and the resulting charge distribution 
depends on the details of pi. For instance, if one sets 
p(q, q') to 6(q)8(q') the probabilities Hf(q) do not depend 



E/(g) =^S(q-N) [ d<pP(N, </>) . (18) 

N J 

A similar result for a simple Josephson junction model 
was cited in . A more general choice of pi preserves the 



^-dependence. One can summarize the situation by say- 
ing that N and 4> are related to canonically conjugated 
variables q and x, that hampers their simultaneous mea- 
surement. 

Since p{x,p) are not positive in general, the P(N, </>) do 
not have to be positive. It might seem that these "nega- 
tive probabilities" can not be measured. Fortunately, it is 
not so. To understand this, let us see how one would mea- 
sure P(N) in the classical case. The only exact way is to 
make use of Eq. (p"7|). One thus measures IT / separately 
and then obtains P(N) from a deconvolution procedure: 
the Fourier transform of P is the ratio of Fourier trans- 
forms of II's. Our main result is that the same decon- 
volution procedure can be applied to Eq. (|l6|), resulting 
in 

P(N d>) = f ^X e m x Pf(<l> + x/2^-x/2) (1Q) 
niV ' 0j J 2n 6 Pl (0 + X /2,0-x/2) • 11 

Since off-diagonal entries of the density matrix can not 
be measured, this expression is not directly applicable. 
In (lj] a scheme was proposed, how this can be circum- 
vented by a repeated measurement of differently prepared 
initial density matrices. This allows one to characterize 
and measure P(N, </>), whatever sign it has. 

In conclusion we have studied the statistical proper- 
ties of supercurrent in short constrictions. An exten- 
sion of the Keldysh technique to account for full count- 
ing statistics of systems containing superconductors was 
developed. In the case of the supercurrent through a 
short constriction (point contact or tunnel junction) the 
cumulant generating function can be found quite gener- 
ally. It shows that charge transfer occurs in units of e 
and 2e, which largely enhanced probabilities in the case 
of contacts with large transmission. The charge trans- 
fer occurs in long trains of electrons passing through the 
contact in either direction, which is a signature of the 
coherent nature of the supercurrent. The relative proba- 
bility of trains in the two directions is determined by the 
thermal occupation and the switching rate between them 
by the broadening parameter in the bulk of the supercon- 
ductors. For tunnel junctions or point contacts between 
gapless superconductors we find the occurrence of nega- 
tive values of P(N, 4>), which questions the interpretation 
as a probability. Accounting for the full time evolution 
of an (idealized) measuring device, we have shown that 
these negative values can indeed be observed. 
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